Study of the "^He crystal surface 
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The evolution of the meniscus of a helium crystal near the (0001) face is traced during a change 
in the boundary conditions at the chamber wall in the temperature range 0.5-0.9 K. The critical 
behavior of the contact angle is studied. An anisotropy is detected in the crystal-glass interface 
energy. New data on the temperature dependence of the elementary-step energy are obtained. 
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In this work, we continue to study the behavior of the 
contact angle that appears when the interface of two con- 
densed phases of '^He (crystal-superfluid) reaches a solid 
wall. The formulation of the problem and a detailed de- 
scription of our optical technique were given in our earlier 
works \m, S] ■ Recall only that the technique consists in 
photographing a crystal using parallel light followed by 
digital processing of images. To improve the images, we 
substantially modified the cryostat design. In the mod- 
ified version, an optical tract passes through windows 
located in the vacuum space of the cryostat and does not 
meet liquid helium. As a result, we were able to decrease 
the image noise by several times. 

Figure[l] schematically shows the cross section of the 
experimental chamber with a crystal at the bottom. Re- 
call also that, when analyzing the meniscus profile near 
the chamber walls, we consider a two-dimensional prob- 
lem in the section plane, since the longitudinal cham- 
ber size (29 mm) is significantly larger than the capil- 
lary constant A ~ 1 mm. The crystal is oriented so that 
the (0001) basal plane is horizontal. The x, z plane of 
a Cartesian coordinate system coincides with the sec- 
tion plane and axis z is vertical. Due to the chosen 
configuration of the side walls, the interface has an 5- 
like shape described by a Z{x) function. Then, angle 
9 = arctg(dZ / dx) determines the surface inclination at 
a given point with respect to the (0001) plane. At the 




point of the end of surface at the wall, we have 9 — Oq, 
and ^0 is one of the quantities to be measured in our 
experiment. The second quantity to be measured is the 
angle of wall inclination tp to axis Ce of the crystal at the 
contact point. 

Angle of contact with the right wall flfl is 



^B. = Tp — ' 
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It is related to the surface energy of the crystal a(0), 
crystal- wall energy Eg, and liquid- wall energy ei by the 
equation 



a{9m)coiiVLii + a 
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where Ae = Eg — £i 

The angle of contact with the left wall is 



^L = iP + 9lo + ij- 
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FIG. 1: Cross section of the experimental chamber 



The boundary condition near the left wall has the form 
a{9Lo)cosnL -a'g{9Lo)smnL = -Ae. (4) 

The right and left walls are made of the same mate- 
rial (polished glass). Therefore, the same physical state 
should take place at both walls at the same values of ip 
for the horizontal orientation of the basal plane in the 
crystal due to the presence of screw twofold axis C2 in 
the symmetry group of the crystal. When C2 rotates, an- 
gle 9lo changes its sign and we will present experimental 
data for the general function 9o{ip) = 9iio{'4)) = — 6'lo(V')- 

To control the possibilities of our optical technique, 
we used the results of digital processing of photographs 
taken from two test samples. Figure[5] shows shape Z{x) 
of the meniscus of liquid"* He filling the lower part of the 
chamber. Here, we also present the surface shape calcu- 
lated in the low-angle approximation 6' ^ 1, 



Z*{x) = A + B cosh-, A 

A' 



where a; = 0,28 erg/cm^ — is the surface tension of liq- 
uid "He atr == 2,3 K [I, pi = 0,145 g/cm^ — is its den- 
sity, and g — is the gravitational acceleration. Fig- 
ure[2]3 shows the dependence SZ{x) = Z{x) ~ Z*{x), 
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FIG. 2: a) Shape Z{x) (line) of the meniscus of liquid ''He and 
(circles) calculated profile, b) Deviation of the experimental 
from the calculated curve. 
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FIG. 3: a) Shape Z{x) of a fiat growing face of a 'He crystal, 
and the dashed line shows the face direction, b) Deviation of 
the experimental curve from a straight line. 



which characterizes the deviation of an experimental 
curve from the calculated curve corresponding to the op- 
timum choice of adjustable parameters A audi?. 

Similar data are presented in Fig. [3] for the crystal the 
(0001) plane of which is close to the horizontal. This 
photograph was taken at a temperature of0.88K, i.e., 
well below roughening temperature Tr of the basal plane. 
Under these conditions, the interface coincides with the 
growing face. In Fig. [3^, dependence Z{x) is approxi- 
mated by a straight line Z* = A + Bx. The difference be- 
tween the experimental results and the optimum straight 
line is shown in Fig.|3jD. In both cases (Figs.I^bilSb)! the 
deviations from the real shape are random (about ±2 /im) 
and are the main origin of errors in determining the sur- 
face shape. The image of the liquid helium meniscus 
can be used to detect the horizontal direction accurate 
to ^ 2 • 10^^ rad. Using photographs of a growing face, 
its direction in the xz plane was determined accurate 
to- 2- IQ-'^rad. 

At the initial stage of experiments, we grew a small 
— 1 mm'' crystal whose basal plane was close to the hor- 



izontal. This procedure was described in detail in [i|, 0] . 
The crystal was then grown slowly to reach the level of 
the glass walls. The growth rate did not exceed — 1 /im/s 
and was stabilized by controlled heating of the ballast 
volume outside the cryostat. The pressure in the ballast 
volume exceeded the equilibrium pressure by '^ 0.4mbar. 

After the chosen temperature was stabilized, we took 
photographs of the crystal profiles at various levels. The 
transitions to the next levels were performed by sequen- 
tial surface melting or growing of the crystal. The differ- 
ence between neighboring levels was 0.2-0.5 mm. 

Once the next level was reached, the chamber with the 
crystal was disconnected from the ballast volume with 
a valve located outside the cryostat at room tempera- 
ture (unfortunately, the design of our device does not 
imply the existence of a "cold" valve). As a result, quasi- 
equilibrium conditions, which are accompanied by a very 
slow helium flow from the filling line to the chamber be- 
cause of a low level of liquid helium in the cryostat bath 
at 4.2 K, were established. Under these conditions, the 
crystal grew slowly at a rate of '-^ 5 ■ 10~^ /xm/s and had 
almost the same shape. 

We now present the measurement results for two crys- 
tals. In the first sample, the basal plane was inclined at 
an angle of 2 • 10~^ rad in the transverse direction (with 




FIG. 4: Series of crystal menisci at T = 0,61 K. The face direc- 
tion is indicated by the heavy dashed line. The photographs 
were taken during melting. 



respect to axis x) and at an angle of 1, 2 • 10~^ rad in the 
longitudinal direction (with respect to axis y). For the 
second sample, the transverse inclination was 8 ■ 10""' rad 
and the longitudinal inclination was 6 • 10"'* rad. After 
measuring the longitudinal face inclination, the optical 
bench was inclined at the measured angle; as a result, 
the error in parallelism between the optical axis and the 
crystal face was corrected. 

The first sample was photographed during gradual sur- 
face melting, and the second sample was grown gradually. 
In both cases, photographs were taken in 20 and40min 
after the chamber was closed with the valve in order to 
control surface relaxation. For the first sample, we took 
two scries of photographs corresponding to temperatures 
of 0.89 and 0.61 K; for the second sample, three series 
of photographs were taken at temperatures of 0.9,0.72, 
and 0.53 K. The results of processing two of the five se- 
ries are shown in Fig. 2] (first sample, T = 0,61 K) and 
Fig. [5] (second sample, T ~ 0,9 K). The scale in the ordi- 
nate axis is higher than that of the abscissa scale by more 
than 60 times. All curves arc located as close as possible 
to each other without conserving the vertical scale. The 
real vertical distance between the centers of the upper 
and lower profiles is about 5 mm. 

The situation seems to be ambiguous from the stand- 




point of an equilibrium surface shape. 

An equilibrium surface shape is known to correspond 
to the minimum of the surface and gravitational energy, 
and the surface rigidity plays a key role in this case. In 
our case, we are dealing with so-called longitudinal rigid- 



ity a 



a" . The surface rigidity was studied by vari- 



ous methods [^tOI • The result substantial for us consist in 
the following. In the temperature range 0.4 K <T < Tji, 
the surface rigidity has a weak anisotropy for atomi- 
cally rough surface regions at sufficiently low angles of 
inclination (0< 0,1 rad). The rigidity is temperature- 
independent {ctQ « 0,245 erg/cm^). A sharp decrease in 
a was only detected [3] at temperatures below 0,3 K and 
low angles 9 < 0,04 rad. 




FIG. 5: Series of crystal menisci at T — 0,9 K. The face direc- 
tion is indicated by the heavy dashed line. The photographs 
were taken during growth. 



FIG. 6: Profile 16 (line) at T = 0,61 K and (circles) approxi- 
mation of the meniscus. 

For all series of our experimental data, surface regions 
at not very low angles of inclination (10~^ ^ ^ ^ Oil rad) 
obey the equilibrium equation aZ'^^ = A.pZ, where 
a = 0,245 erg/cm^, and Ap is the difference between the 
densities of the solid and liquid phases. Fig. 6 shows 
the approxiination of the profile 16 by the solution to 
Eq. (??), Z = ^sinh^;^, where A = {a/ApgY^^. The 
dashed line shows the measured direction of the (0001) 
face inclined at 2 ■ 10~^ rad to the horizontal. 

Note that the meniscus shape can also be de- 
scribed [2| with a satisfactory accuracy by the function 
Z = ±(x — xo)^b~'^, where 6 ~ 1 mm, which corresponds 
to the standard theory of a vicinal surface as the ech- 
elon of steps that repel each other as oc x~^. However, 
this picture is in conflict with the results of measuring 
a from the spectrum of crystallization waves at not very 
low angles. 

Similar results of processing profllc 32 Fig. [5] are 
shown iuFigl?] In this case, the left and right seg- 
ments of the profile are approximated by the formula 
Z± = A sinh — -r-^ at the same amplitude A and param- 
eters x± , differing by « 1 mm. 

Thus, we can state that the surface is in equilibrium 
at not very low angles 10"^ < 6 < 10~^. The situation is 
radically different at angles 6 < 10~^ rad: non-monotonic 
profiles often appear, which indicates the absence of equi- 
librium (sec Figs. l5|4|) . The lengths of such regions along 
axis X change within several millimeters, and the irreg- 
ular deviations from the vertical are 10 20 /xm. Such 

phenomena are usually related to lattice defects, mainly 
dislocations, and we cannot exclude this explanation. We 
only recall that the crystals were grown with all precau- 
tions necessary in such cases. 
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FIG. 7: Profile 32 (line) at T = 0,9 K and its approximation 
(circles). 



However, the shapes of some surface fragments can 
hardly be explained by the presence of dislocations. For 
example, the center portion of profile 23 (Fig. [5]) con- 
tains a linear (within the limits of experimental error) 
^ 4 mm segment inclined at an angle of 0.003 rad to the 
(0001) face. The next frame (profile 24) contains a flat 
region reaching the left wall and inclined at an angle of 
0.01 rad to the opposite side. The presence of such ex- 
tended surface regions is thought to be hardly explained 
by the existence of defects in the crystal volume. Unfor- 
tunately, we have no certain considerations regarding the 
nature of these metastable states. It is difficult to draw 
any quantitative conclusions concerning the equilibrium 
meniscus shape (and the contact angle) and, hence, the 
angular dependence of the surface rigidity at low angles 
under these conditions. 

The authors of Q theoretically predicted a phe- 
nomenon caused by the jump of derivative a^ at 6 — 0. 
This phenomenon consists in the fact that the state where 
an atomically smooth face is in immediate contact with 
the wall takes place in the angular range '0_ < -0 < -0+, 
determined by the relation 



Ae — ao sin ip\ < [3 cos i/j 



6 = +0. In this case, a plateau 
o(V-') dependence in a certain an- 



where (3 — a'g > at 
should appear in the 
gular range. Figure|S] shows the results of processing our 
photographs for this dependence. The plateau is seen to 
exist. However, the Oo{ip) dependence differs substan- 
tially from the root behavior suggested in [3]. The con- 
tact angle approaches the plateau linearly at not very low 
angles Oq ^ 1°. 

Recall that the meniscus profile was determined from 
an analysis of a diffraction pattern, as was described in 
our earlier works [ij, i2i]. However, it is difficult to per- 
form this analysis near the crystal-liquid-wall contact 
line, i.e., in the region where the contact angle is to 
be determined. Here, the diffraction pattern is compli- 
cated because of the approaching of two relatively simple 
diffraction patterns from the wall and the crystal-liquid 




FIG. 8: Dependence 6o{'tp) for two series of profiles obtained 
at a; r = 0,61 K b) T = 0,9 K. Open and solid symbols cor- 
respond to the right and left sides of a meniscus, respectively, 
(solid line) Calculation by Eqs. Q,®,®, and (O. 



interface at a distance of w 0.3 mm, which is compara- 
ble with A w 1.2 mm. The value of Oq at 9q > 0,01 rad 
under these conditions was determined by the extrap- 
olation of the derivative of the function approximating 
the meniscus Z = A sinh ^~^° to the calculated point of 
contact with the wall. The region of low contact angles 
6*0 < 0.01 rad, where the root law is likely to be valid, 
requires an additional investigation. 

The observed behavior of the meniscus at angles 
10^^ < < 10^^ rad can be explained under the follow- 
ing assumptions: function a{0) at these angles has the 
form 



ao + P\e\+a'^ 



2 ' 



(5) 



and solid helium-glass interface energy Ss is a func- 
tion of angle ■0. Moreover, we believe that, in the 
temperature range under study, only the step energy 
changes with temperature and the other parameters 
ao « 0,172 erg/cm^ Q, «(,' = ao - ao « 0,073 erg/cm^, 
and function As{tp) almost reached their values charac- 
teristic of zero temperature. 

Figurc[9] shows the Ae('0) function calculated by Eqs. 
(21) and dH) using the experimental data for the two series 
of measurements. Parameters /3 for each series were cho- 
sen so that the imaginary extensions of the Ae(i/') func- 
tions calculated at -0 < 'ip_ and ^p > V'+, arc matched in 
the plateau region. The results for the Ae(-0) functions 
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FIG. 9: Dependence Ae{tp) plotted from the data of the 
two series of measurements at T = 0,9 K (open symbols) and 
(solid symbols) T = 0,61 K (solid line) Calculation by Eq. ©. 
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FIG. 10: Dependence /3(r). 



thus obtained coincide for both series within the hmits of 
experimental error. We take into account the symmetry 



of the helium crystal, neglect the azimuthal anisotropy, 
assume that Ae{ip) is an analytical function of angles, 
and parametrize it as (in erg/cm^) 

Ae(V') = 0,128 -0,013cos27/' + 0,022 cos 4^^, (6) 

The calculated 9o{'ip) dependence with allowance for the 
anisotropy of Ae(^) (see Eq. ([6])) are shown as the solid 
lines in Fig 151 

Using this scheme and Eq. dH), we then determined the 
values of (3 for the other series. Note that the neglect of 
anisotropy £s('0) leads to an increase in the estimate of /3 
by about 30%. The obtained temperature dependence of 
/3 is shown in Fig.^Dl (solid circles) . Here, we also present 
the results from [7[ (open circles), which were obtained 
in the temperature range 50—250 K by an analysis of a 
spectrum of crystallization waves, and from [J| (open tri- 
angles), which were obtained by an analysis of the (0001) 
face growth kinetics at temperatures close to the rough- 
ening temperature, and the value of f3 obtained in [ij 
(cross) at a temperature of 0.72 K. 

Thus, our data on (3 characterizing the surface energy 
at angles 0,01 < < 0, Irad agree with the data in 7| 
obtained at 0<O, Olrad. However, the authors of [7| 
detected a sharp decrease in a at temperatures lower than 
0.25 K and angles 9 < 0,01rad, which agrees with the 
theoretical concepts of vicinal surfaces, and we detected 
non-analytical contribution f3\9\ to the surface energy in 
the angular range 0,01 < 9 <0, Irad, which cannot be 
explained theoretically at finite and weakly temperature- 
dependent a. 
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